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Abstract 



The origin of the bulk of cosmic rays (CRs) observed at Earth is the topic of a century long investigation, paved with successes 
and failures. From the energetic point of view, supernova remnants (SNRs) remain the most plausible sources of CRs up to rigidity 
~ 10 6 - 10 7 GV. This confidence somehow resulted in the construction of a paradigm, the so-called SNR paradigm: CRs are 
accelerated through diffusive shock acceleration in SNRs and propagate diffusively in the Galaxy in an energy dependent way. 
Qualitative confirmation of the SNR acceleration scenario has recently been provided by gamma ray and X-ray observations. 
Diffusive propagation in the Galaxy is probed observationally through measurement of the secondary to primary nuclei flux ratios 
(such as B/C). There are however some weak points in the paradigm, which suggest that we are probably missing some physical 
ingredients in our models. The theory of diffusive shock acceleration at SNR shocks predicts spectra of accelerated particles 
which are systematically too hard compared with the ones inferred from gamma ray observations. Moreover, hard injection spectra 
indirectly imply a steep energy dependence of the diffusion coefficient in the Galaxy, which in turn leads to anisotropy larger than 
the observed one. Moreover recent measurements of the flux of nuclei suggest that the spectra have a break at rigidity ~ 200 
GV, which does not sit well with the common wisdom in acceleration and propagation. In this paper I will review these new 
developments and suggest some possible implications. 
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1. Introduction 



aspects of the theory of the origin of CRs that I think 



Understanding the origin of cosmic rays means to be 
able to unfold a complex chain of physical processes 
that go from the acceleration of a few charged parti- 
cles in potentially unknown sources to their propagation 
in the interstellar medium (ISM). Each of these steps 
involves plasma physics processes often in a non lin- 
ear regime. The whole picture is usually made fuzzier 
by the poor knowledge of the values of environmen- 
tal parameters (e.g. densities, temperatures, fraction 
of ionized material). Moreover most observations in- 
volve quantities (such as fluxes, chemical composition, 
anisotropy) measured at Earth that are averages over 
long (propagation) times and over numerous sources, 
potentially very different. It is not surprising that a cen- 
tury after the discovery of cosmic rays we are still de- 
bating about some aspects of this problem. The purpose 
of this short review is to provide a biased view of some 



are understood and some that are not understood or are 
poorly understood. I will limit myself to CRs that we 
think originate inside the Galaxy, though the definition 
of this transition energy is somehow part of the aspects 
that need to be discussed. 

So far, the only theoretical framework that reached 
a sufficient level of elaboration to deserve the name 
"model" is the one based on SNRs as the main sources 
of the bulk of Galactic CRs. In fact, the lack of a rea- 
sonable alternative has elevated the model to the rank of 
a paradigm. CR acceleration is believed to occur at the 
forward shock of SNRs through diffusive shock accel- 
eration (DSA). The test particle theory of DSA (see fH 
for a review) predicts that the spectrum of accelerated 
particles is N(E) cc E~ y with y — (r + 2)/(r - 1), where 
r is the compression factor of the shock. For typical 
parameters of SNR shocks, r»4 and y — 2. The prop- 
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agation of CRs in the Galaxy is usually parametrized 
through a diffusion coefficient D(E) oc E s . For nuclei 
for which spallation is negligible the equilibrium spec- 
trum observed at Earth is n(E) oc N(E)/T esc (E), where 
T esc {E) ~ H 2 /D(E) ~ E~ s , where H is the size of the 
halo of the Galaxy. It follows that n(E) oc E~ 7 ~ 5 . Com- 
parison with observations leads to 6 ~ 0.7 if y » 2. 
6 ~ 0.6 - 0.7 is also compatible with the low energy 
observed slope of the B/C ratio. This might look like 
a self-consistent picture at first sight, but two problems 
immediately arise: 1) the observed gamma ray spectra 
in the 1-100 GeV range suggest a steeper CR spec- 
trum, with y ~ 2.3 - 2.4 [2j; 2) A diffusion coefficient 
D(E) oc E would lead to a large scale CR anisotropy 
much larger than observed 0. The non-linear theory 
of DSA (NLDSA, see |@) for a review) makes these 
two problems even more severe: the theory accounts for 
the dynamical reaction of accelerated particles on the 
shock, which is responsible for the formation of a pre- 
cursor upstream of the shock. In NLDSA the compres- 
sion factor felt by accelerated particles becomes a func- 
tion of energy, which reflects in concave spectra, which 
turn out to be even harder than E~ 2 at E > 10 - 100 
GeV. As shown in [5 1, injection of CRs in SNRs through 
NLDSA leads to require D(E) oc E 0J5 , that is hardly 
compatible with the observed large scale anisotropy. 

It is clear that the problems of acceleration in the 
sources and propagation in the ISM are tightly con- 
nected with each other and cannot be studied as two in- 
dependent problems. For instance, the assumption that 
the diffusion coefficient D(E) is a given function of en- 
ergy is most likely inappropriate, in that CRs themselves 
are able to create the scattering centers responsible for 
their diffusive motion. In turn, the diffusion coefficient 
affects the CR spectrum observed at Earth. This phe- 
nomenon, discussed by (6] [7j might play an important 
role in understanding features on the CR spectrum re- 
cently reported by PAMELA QD and CREAM g). 

This paper is organized as follows: in ^2] I will dis- 
cuss the implications of the NLDSA in terms of mag- 
netic field amplification and how this might mitigate the 
spectral problem illustrated above. In ^3] I will discuss 
the important role of CRs to self-generating the scatter- 
ing centers on which diffusion in the ISM takes place. 
The implications for nuclear spectra will also be illus- 
trated. In ^4] I will discuss an important recent devel- 
opment in the theory of CR acceleration in SNRs, in 
connection with the presence of neutrals hydrogen in 
the acceleration region. In ^5]l will comment on the im- 
plications of the SNR paradigm for the transition from 
Galactic to extragalactic CR. A summary will be pre- 
sented in §j6] 



2. Magnetic field amplification, Maximum energy 
and velocity of scattering centers 

NLDSA describes the process of particle acceleration 
at collisionless shocks relaxing the test particle assump- 
tion, namely taking into account the CR dynamical re- 
action on the shock. Moreoever, J3] [10] QT) also in- 
troduced in the theory the phenomenon of CR-induced 
magnetic field amplification and the dynamical reaction 
of the magnetic field on the background plasma. 

The magnetic field can be amplified because of the 
streaming instability induced by the super-alfvenic mo- 
tion of CRs in the upstream plasma Ifl2l [131 [T4l . The 
instability leads to resonant growth of modes with 
wavenumber k ~ l/ri, where ri is the Larmor radius of 
particles that dominate the particle number at the shock. 
The waves produced by accelerated particles can also 
resonantly be absorbed by the same particles thereby 
leading to their diffusive motion. More recently ifTBI 
found that non resonant modes may grow much faster 
than the resonant modes on scales k » Although 
this phenomenon may be relevant for magnetic field am- 
plification, it is probably not the main mechanism for 
scattering the particles and for their acceleration at the 
shock. 

The amplification of magnetic field due to the stream- 
ing of CRs with the shock was first introduced in the 
theory of NLDSA by ||T7l[l8|. Its effect for achieving 
higher maximum energies was later investigated in [ 19 1. 
Magnetic field amplification is also needed to explain 
the narrow X-ray rims detected in virtually all young 
SNRs (see |20| and references therein for a recent re- 
view). 

Finally the creation of turbulent magnetic field in the 
shock region might have an important effect in deter- 
mining the shape of the spectrum. As first emphasized 
by IfTZi l. the relevant shock compression factor for par- 
ticle acceleration is the ratio of velocity of the scatter- 
ing centers. If the waves are slow in the frame comov- 
ing with the plasma, then this is basically the same as 
the compression factor of plasma speeds. On the other 
hand, in the non-linear regime, the waves might acquire 
a substantially higher speed (perhaps of the order of the 
Alfven speed calculated in the amplified field). In this 
case the relevant compression factor may be somewhat 
smaller and reflect into steeper spectra of accelerated 
particles Il23l l24l . It is worth recalling that the details 
of this phenomenon are very model dependent: whether 
the spectrum becomes harder or softer depends on the 
helicity of waves, which is hard to predict, especially 
in a non-linear regime such as the one that is expected 
close the shock front. In |f2l [24) the authors assume 
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Figure 1 . Spatially integrated spectral energy distribution of the Tycho 
SNR 1 29 1 - The curves show the calculated multifrequency spectrum 
from the radio to the gamma ray band. 



that wave velocity can be estimated as the Alfven speed 
as calculated in the amplified field, vw - SB/ 
with SB » Bq, Bq being the pre-existing magnetic field. 
The situation might be appreciably more complex than 
that: as discussed in lfl4l [131 [TBI the growth rate of 
the resonant modes excited by CRs in the regime of 
high CR acceleration efficiency have a phase velocity 

V <P ~ ( ! « ZvjC ) ' where uqr is the density of accelerated 
particles, n ; is the ion density, v s the shock velocity and 
c the speed of light. It is easy to see that may easily 
exceed the Alfven speed. In this regime the growth rate 
of the waves also changes. 

The simple recipe of ||23l l24l was used to calculate 
the spectrum of CRs at the Earth and it was shown that 
the required Galactic diffusion coefficient is D(E) oc 
E 0M (see also [26 1), which alleviates but does not solve 
the anisotropy problem. The finite velocity of the scat- 
tering centers was also used to calculate the multifre- 
quency spectrum of the Tycho SNR |29|. The results of 
this calculation are illustrated in Fig. [T] where I show the 
multifrequency spectrum from the radio band to gamma 
rays. The relatively steep spectrum in the gamma ray 
range, that provides a good fit to the data points from 
Fermi-LAT Il27ll and Veritas l28l , illustrates well the ef- 
fect of the velocity of scattering centers (the gamma ray 
spectrum vF v would be roughly flat in the absence of 
this effect). The predicted strength of the magnetic field 
reproduces the brightness profile of the non-thermal X- 
ray emission, as well as the synchrotron spectrum from 
radio to X-rays. The case of Tycho represents the first 
convincing instance of a SNR accelerating CRs to ener- 
gies of the order to > 500 TeV. 

It is worth noticing that the magnetic field amplifi- 



cation as observed in the form of narrow X-ray rims in 
SNRs might also be the result of purely hydrodynam- 
ical processes [30 1, as due to the presence of density 
inhomogeneities upstream that result in vorticity at the 
shock crossing and magnetic field amplification due to 
the wrapping of the eddies downstream. This mecha- 
nism operates downstream of the shock and no mag- 
netic field amplification should be expected upstream. 
This is a very important point in that DSA requires ef- 
fective scattering of the particles on both sides of the 
shock and in fact if magnetic field amplification only 
takes place downstream, no significant increase in the 
maximum energy should be expected. One noticeable 
exception to this statement appears if the shock is quasi- 
perpendicular in which case particle acceleration may 
proceed through drifts and may in principle be very fast. 
In this perspective the case of Tycho is especially impor- 
tant, since it is a supernova of Type la that exploded in 
the normal ISM and preferential acceleration at perpen- 
dicular shocks should result in a bilateral structure that 
is quite unlike the spherical appearance of the Tycho 
SNR. 

3. Propagation of CRs in the Galaxy 

The propagation of CRs in the Galaxy is best un- 
derstood in terms of diffusion in a disordered magnetic 
field, as shown by the relative abundance of secondary 
nuclei and their parent primary nuclei (most notably 
boron versus carbon), which indicate confinement times 
in the Galaxy volume which exceed any ballistic time 
by orders of magnitude. Similar conclusions can be 
reached based on the abundances of some unstable iso- 
topes (such as m Be). The propagation of nuclei and 
leptons in the Galaxy is usually modeled with the help 
of numerical codes such as GALPROP and DRAGON. 
Although a substantial agreement is claimed between 
predictions and data, usually the price to pay in these 
models is that of imposing breaks in the diffusion coef- 
ficient and/or the injection spectrum in order to repro- 
duce observations. Moreover there is some degeneracy 
between different models such as models with reaccel- 
eration and D(E) oc E lli and models without reaccel- 
eration and D(E) oc E 05 . The physical origin of these 
alleged breaks is unknown, which makes the fit of the 
data somewhat unsatisfactory from the physical point of 
view. The recent data by PAMELA and CREAM have 
made the situation more puzzling. 

The PAMELA [ 8 ] experiment has provided evidence 
that the spectra of CR protons and helium nuclei have 
a break at rigidity ~ 200 GV. The spectrum of protons 
measured by PAMELA has a slope 2.85 + 0.Q\5{stat) ± 
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0.004(sysf) at R < 240 GV and 2.67 + 0.03 + 0.05 at 
R > 240 GV. The spectrum of He nuclei has slopes 
2.766 ± 0.01 ± 0.027 and 2.477 + 0.06 + 0.03 in the 
same rigidity ranges. The CREAM experiment mea- 
sured these spectra at higher energies and found slopes 
compatible with those of PAMELA for/? > 200 GV (for 
protons the high energy slope is 2.66 + 0.02 and for He 
the slope is 2.58+0.02). This finding suggests that some 
new physical phenomenon might be occurring at a scale 
of ~ 200 GV either in the physics of acceleration or in 
the diffusion of CRs throughout the ISM. A few models 
have been proposed to explain the spectral hardening: 
OTI suggested that protons and He may be injected with 
slightly different spectra as a result of a different scaling 
of their injection rates during the temporal evolution of 
the SNR. This idea is very important to understand why 
the two spectra are slightly different but does not ad- 
dress the issue of the spectral hardening observed for 
both p and He nuclei at ~ 200 GV. In E3, the author 
points out that a spectral hardening may be understood 
if the diffusion coefficient in the halo is different from 
that in the disk. In |7] the spectral hardening was at- 
tributed to the interplay between the self-generated tur- 
bulence induced by CR streaming in the ISM and the 
non linear Landau damping that induces the cascade of 
turbulence from a large scale injection, possibly due to 
SN explosions. The problem of diffusion of CRs in this 
partially self-generated turbulence is non linear and can 
be solved only by coupling the transport equation with 
the equation describing magnetic field generation and 
damping. A transition in the diffusion properties of CRs 
is predicted by [7] to appear at rigidity ~ 200 GV: the 
diffusion coefficient at low energy has an approximate 
scaling D(R) ~ E 0J , while at R > 200 GV, the scaling 
becomes D(R) ~ E 1 ^ for a Kolmogorov phenomenol- 
ogy (see Fig. [2}. At energies below ~ 10 GeV prop- 
agation becomes dominated by advection of CRs with 
waves moving at Alfven speed away from the disc of the 
Galaxy. Since the self-generated waves are produced 
along the CR gradient, there are no waves moving to- 
wards the disc, hence one should not expect appreciable 
reacceleration in this model. The calculation of [7| was 
also used to show that this approach can accomodate the 
spectrum of CRs inferred by 1 33 34 1 using the gamma 
ray spectrum observed by Fermi-LAT from the direction 
of clouds in the Gould belt. The calculations of prop- 
agation of other nuclei including the effects of spalla- 
tion of nuclei and wave generation and absorption from 
all nuclei is currently being done (Aloisio and Blasi, in 
preparation). 

The main lesson to be learnt from this model is that 
the ability of CRs to generate the conditions for their 
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Figure 2. Diffusion coefcient induced by streaming instability of CRs 
and cascading from a large spatial scale of 50 pc (figure from Ref. 
03). 



diffusion is a widespread phenomenon, and it is likely 
to be very important in the acceleration sites as well as 
in the propagation of CRs in the ISM. 

4. DSA in the presence of neutrals 

It is important to realize that SNR shocks are colli- 
sionless, namely their formation is not due to particle- 
particle scattering but rather to the mediation of elec- 
tromagnetic instabilities. The thickness of the shock 
front is expected to be of the order of the Larmor ra- 
dius of thermal protons behind the shock. Interest- 
ingly, even electrons and protons, due to their different 
mass are expected to thermalize to different tempera- 
tures, T e /T p ~ m e /m p , behind the shock. Other colli- 
sional and collisionless processes may equilibrate elec- 
trons and protons downstream. Neutral atoms that cross 
a collisionless shock do not experience any jump but 
are coupled to the background plasma through the pro- 
cesses of charge exchange and ionization Il35ll36ll . 

The presence of partially ionized material in the ac- 
celeration region may profoundly change the way DSA 
works, mainly because of 1) ion-neutral damping of 
waves which may stop the growth of CRs induced 
waves and hamper the acceleration process, and because 
of 2) the dynamical reaction of neutral material in prox- 
imity of a collisionless shock front. In addition, the 
presence of neutral atoms may provide us with a pre- 
cious diagnostic tool of the acceleration process, as we 
discuss below. 

The dynamical reaction of neutrals on the shock is 
mainly due to the phenomenon of the neutral return flux 
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Figure 3. Slope of the spectrum of accelerated test particles for E = 1 , 
10, 100, 1000 GeV, as a function of the shock speed (figure from Ref. 
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(NRF) first investigated in the context of SNRs by 1571 . 
The main coupling between neutrals and ions at a col- 
lisionless shock is due to charge exchange and ioniza- 
tion, that are activated when the net relative speed be- 
tween the two components is non-zero. Downstream of 
a collisionless shock ions are slowed down and heated 
up, while neutrals cross the shock and keep their ini- 
tial velocity. The charge exchange reactions occur- 
ring in this situation may eventually produce neutrals 
moving with large bulk velocity in the direction of the 
shock and these ions may recross the shock toward up- 
stream. Charge exchange and ionization with the up- 
stream plasma lead to deposition of energy and mo- 
mentum of these neutrals upstream, which in turn re- 
sults in heating of this gas. For shock speed < 4000 
km/s this NRF considerably changes the structure of the 
shock and leads to the formation of a shock precursor, 
similar to the one induced by CRs but in general on a 
smaller spatial scale. In 0711 the authors showed that the 
spectrum of test particles accelerated at such shock may 
visibly deviate from the standard predictions of DSA 
and account for steeper spectra of accelerated particles. 
For shocks faster than ~ 4000 km/s, on average a neu- 
tral crossing the shock gets ionized before suffering a 
charge exchange, therefore the NRF is suppressed. The 
net effect of the NRF is to heat the plasma in the pre- 
cursor so that the Mach number of the shock decreases, 
namely the compression factor gets smaller. Since the 
spectrum of particles accelerated at the shock depends 
on the compression factor, the NRF leads to steeper par- 
ticle spectra. Clearly the effect is the largest for particles 
diffusing on spatial scales comparable with the scale of 
the neutral induced precursor. Since the diffusion coef- 



ficient is a growing function of energy, the effect of the 
NRF shows more prominently at lower energies. In Fig. 
[3]we show the slope of the spectrum of accelerated test 
particles as a function of the shock velocity for different 
values of the particle energy. The ionization fraction is 
assumed to be 50%. As expected, the spectral steepen- 
ing is mostly present for v S h < 3000 km/s, and it extends 
to particle energies of ~ TeV. At high enough energy 
the standard slope y — 2 is recovered. 

As mentioned above, neutral hydrogen atoms in the 
acceleration region also provide us with an important 
diagnostics of the acceleration process, through their 
Balmer emission. The Balmer line emission from hy- 
drogen in the shock region is becoming a powerful 
tool to measure the CR acceleration efficiency in SNR 
shocks. The idea is relatively simple: it is well known 
that the Balmer line produced by neutrals that suffered 
charge exchange with hot ions downstream has a width 
that reflects the temperature of ions, while neutrals that 
did not suffer charge exchange emit a narrow Balmer 
line with width ~ 20 km/s, corresponding to the up- 
stream T ~ 10 4 K hydrogen temperature. If CR ac- 
celeration is efficient, part of the ram pressure upstream 
of the shock pu 2 is channelled into accelerated parti- 
cles instead of heating. Therefore the plasma tempera- 
ture downstream is lower if particle acceleration is effi- 
cient and the corresponding broad Balmer line becomes 
correspondingly narrower. On the other hand, efficient 
CR acceleration produces a shock precursor upstream 
that leads the ionized plasma to slow down with re- 
spect to neutrals. Therefore charge exchange is also 
activated upstream and the narrow component of the 
Balmer line may become broader. As shown in [38 1, 
for slow enough shocks the NRF also produces a simi- 
lar effect but it results in the formation of a component 
of the Balmer line with width intermediate between the 
broad and the narrow line. 

The appropriate use of Balmer lines as a diagnostic 
tool of CR acceleration must be based on a careful mod- 
eling of the interplay between ions, neutrals and CR. 
A theory of NLDSA in the presence of neutrals is cur- 
rently being developed. A semi -quantitative assessment 
of the CR acceleration efficiency in selected SNRs has 
been made in the past few years: in l39l the authors 
measured the width of the broad Balmer line in a re- 
gion of the SNR RCW86 and found that the inferred ion 
temperature could be interpreted as the result of > 50% 
of the ram pressure being converted to accelerated par- 
ticles. A review of Balmer lines as diagnostics of CR 
acceleration can be found in |40|. 
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5. Implications of the SNR paradigm for the transi- 
tion to extragalactic CR 

The amplification of magnetic field in proximity of 
a SNR shock is crucial to accelerate CRs up to rigidity 
R ~ 10 6 — 10 7 GV f!9\ . In fact higher energies are very 
hard to achieve unless some very special conditions are 
realized (fast quasi perpendicular shocks), as speculated 
in ll24l in connection with SN type Ic. In the minimal 
version of the SNR paradigm, the all-particle spectrum 
of CRs can be explained in terms of particle acceleration 
in core collapse supernovae and Type la supernovae: the 
knee results from a change in chemical composition and 
the Galactic CR spectrum ends with an iron dominated 
composition at E ~ 2 x 10 17 eV. At this point the transi- 
tion to extragalactic CRs should be realized. There are 
basically two models of extragalactic CRs in which the 
transition occurs in this energy region. 

For long time it has been taken for granted that the 
ankle in the CR spectrum, at ~ 10 19 eV, is the spectral 
signature of the transition from a steep Galactic spec- 
trum to a flatter extragalactic spectrum. The situation 
has however changed and the nature of the ankle ques- 
tioned as a consequence of two developments: a) in Ref. 
fiTTl the authors noticed that Bethe-Heitler pair produc- 
tion leaves a distinct feature in the spectrum of CRs 
propagating on cosmological scales. The feature takes 
the form of a dip whose shape fits very well the observed 
modification factor for all experiments, with the possi- 
ble exception of the one measured by the Pierre Auger 
Observatory. In this model CRs with energy > 1 EeV 
are of extragalactic origin, and the transition occurs at 
the second knee, b) In Ref. |42| the authors discussed 
the possibility that UHECRs may be nuclei with a mixed 
chemical composition. In this model the Galactic com- 
ponent of CRs ends at energy ~ 2 EeV. The so-called 
disappointing model introduced in Ref. 11431 is a spe- 
cial case of the mixed composition model, in which the 
maximum energy of protons is relatively low, ~ 4x 10 18 
eV, and the iron spectrum extends to ~ 10 20 eV. The 
model is disappointing in that the flux suppression at 
~ 10 20 eV is not the GZK feature but rather the intrin- 
sic cutoff in the source spectrum and no correlation with 
sources is expected because of the heavy composition at 
the highest energies. 

Both these models (dip and mixed composition) lead 
us to expect that Galactic CRs may end in the < EeV 
region (with a composition dominated by heavy nuclei) 
rather than at the ankle. This conclusion also appears 
to be supported by a recent investigation of cosmic ray 
anisotropies fl4ll . 

The spectrum of cosmic rays is shown in Fig. [4] for 
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Figure 4. Spectrum of CRs in the dip model overlapped to the Galactic 
CR flux. 
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Figure 5. Spectrum of CRs in the disappointing model with the Galac- 
tic CR flux as in Fig. [4] 



the dip model and in Fig. [5] for the mixed composition 
model in its disappointing model configuration. In these 
calculations the spectrum of Galactic CRs was taken 
from P31l . In the right panel, the different lines rep- 
resent the fluxes of different chemicals, the solid line 
representing the total flux. One can clearly see that at 
high energy the chemical composition is dominated by 
heavy nuclei. The transition region is better described 
by the dip model. 

The main discrimination among dip model, mixed 
composition model and ankle model is based upon the 
measurement of the chemical composition, especially in 
the transition region |46|. In the dip model the elonga- 
tion rate is expected to show a sharp transition from a 
heavy Galactic composition at energy below the second 
knee to a light composition at E > 10 18 eV (see |47|). 
The transition is predicted to be concluded at 10 18 eV, 
where the composition is completely dominated by pro- 
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tons. The dip model works well provided the abundance 
of helium in the primaries is smaller than ~ 10% in flux 
of extragalactic CR. 

In the mixed composition model the extragalactic 
CRs are made of a mix of different elements, and the 
flux at Earth is the result of a complex chain of interac- 
tions: photodisintegration of nuclei leads to lighter com- 
position even when starting with pure iron at the source. 
One can change the injection spectra and the compo- 
sition at the source to match the observed composition 
and spectra, although in general it is hard to explain a 
transition to heavy composition at ultra high energies, 
as it is observed by Auger. In the case of the mixed 
composition model the elongation rate shows a gradual 
transition from heavy galactic CRs to somewhat lighter 
extragalactic CRs. The so-called disappointing model 
reproduces the heavy composition of Auger at high en- 
ergy by construction. 

In the traditional ankle model the transition reflects 
in a gradual change from an iron dominated composi- 
tion of Galactic CRs that extend to ~ 10 19 eV to a pure 
proton composition, reached at energies ~ 5 x 10 19 eV, 
so that in this case only a small region of energies is 
filled by the extragalactic CR component. This model 
does not fit the elongation rate as observed in any of the 
current experiments and is not immediately compatible 
with the SNR paradigm for the origin of Galactic CRs. 

6. Summary 

Recent gamma and X-ray observations of SNRs have 
strengthened the case in favor of SNRs as the main 
sources of Galactic CR. However several loose ends re- 
main as they arise when a serious attempt is made to 
make sense of what we measure. On one hand the spec- 
trum of particles generated by DSA for strong shocks 
is very close to E~ 2 , this prediction is at odds with the 
measurement of the large scale anisotropy |3] since it 
requires the diffusion to be exceedingly fast. Moreover 
the spectra of SNRs observed in gamma rays, when in- 
terpreted as a result of production and decay of neutral 
pions, suggest steeper spectra [23]. The difference be- 
tween instantaneous and escaping CR spectrum ll2Tll22l 
does not appear to be sufficient to explain this discrep- 
ancy. The usage of NLDSA instead of test particle 
DSA makes the problem even more severe, although it 
has been proposed that somewhat steeper spectra might 
be produced if the finite velocity of the scattering cen- 
ters is taken into account (23]. This idea has been 
tested versus the multifrequency modeling of the Tycho 
SNR [29] and versus the spectra of propagated nuclei 



in the ISM |26l|24). In the latter case a Galactic diffu- 
sion coefficient D(E) oc E 0M was required (instead of 
D(E) oc E - 15 found by JS | in the absence of this effect). 

An independent support to the SNR paradigm is pro- 
vided by the detection of narrow non-thermal X-ray 
rims in numerous SNRs. The thickness of these fila- 
ments can be best explained in terms of magnetic field 
amplification, possibly induced by the streaming insta- 
bility excited by CRs at SNR shocks. The observed 
fields are of order of 100 - 1000//G and may explain 
the acceleration of CRs up to 10 6 - 10 7 GV rigidities, 
as required by CR data. These findings outline a situa- 
tion in which the Galactic CR spectrum ends at energies 
~ (2 - 3) x 10 eV with an iron dominated chemical 
composition |45|. 

The simple picture that has survived for several 
decades of a simple diffusive propagation of CRs in the 
Galaxy has recently been threatened by some observa- 
tional findings. The PAMELA experiment found that 
both the proton and He spectra show a break at rigidity 
~ 200 GV. The He spectrum is also found to be sys- 
tematically harder than the proton spectrum at all en- 
ergies. The hardening of spectra might be the signa- 
ture of a physical process that affects either accelera- 
tion or propagation in the ISM. In [7 1 the authors sug- 
gest that this change of slope is the result of a tran- 
sition from self-generated waves, where the scattering 
centers responsible for CR diffusion are generated by 
the same CRs through streaming instability, to diffusion 
in a background of waves cascading from larger scales 
for instance through non-linear Landau damping (wave- 
wave coupling). In this model the diffusion coefficient 
is an output of the calculations, and shows a complex 
structure, as illustrated in Fig. [2] The spectrum of CRs 
fits the slopes measured by PAMELA and the transition 
energy between the two regimes of diffusion is naturally 
found to be around ~ 200 GeV (for protons). Interest- 
ingly the slope of the PAMELA spectrum at R < 200 
GV is the same as that found by |[33l [34 1 by analyzing 
the gamma ray spectra of clouds in the Gould belt and 
inferring the parent CR spectrum. The slope of the CR 
spectrum at R > 200 GV is confirmed by measurements 
carried out with the CREAM experiment. 

Although there is no doubt that gamma ray observa- 
tions are playing a crucial role in establishing the role 
of SNRs as the main sources of Galactic CR, it is also 
clear that uncertainties in the environmental conditions 
in which acceleration takes place creates some ambigu- 
ity in the conclusions and often these uncertainties do 
not allow us to discriminate between a leptonic and a 
hadronic origin of the observed gamma radiation. In 
some SNRs different regions may radiate gamma rays 
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due to different mechanisms and the angular resolution 
of the gamma ray telescopes may not be sufficient to 
make these situations clear. Other methods to establish 
whether a SNR is accelerating CRs efficiently should 
be sought after. In this perspective, the most interest- 
ing method is the observation of optical Balmer lines. 
The width of Balmer lines is affected by the presence of 
CRs in the acceleration region: the broad component is 
narrower if an appreciable fraction of the ram pressure 
entering the shock is being channelled into CR acceler- 
ation. On the other hand if CRs are being accelerated 
effectively, a precursor is formed upstream of the shock, 
and this may lead to a broader narrow component of 
the Balmer line. Spatially resolved observations of the 
shape of the Balmer line in SNRs may turn out to pro- 
vide a wealth of new information about the process of 
particle acceleration in SNRs. 
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